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> ■ 

Friedmann-Lemaitre universes driven by a scalar field, spatially closed and bouncing, were recently studied in [1], 
; I ■ with the conclusion that the spectrum of their large scale matter perturbations was generically modified when going 
through the bounce. In this Note we extend this result to a wider class of bouncing scale factors and give the properties 
of the scalar field potentials which drive them. In doing so we throw light on the hypotheses which underlie the models 
and discuss their cogency. 
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I. Introduction 

Since the invention of the string-inspired, "pre-Big Bang" [2] and "ekpyrotic" or "cyclic" [3] universes, there has 
been renewed interest in 4-dimensional, general rclativistic, bouncing, Fricdmann-Lcmaitrc models, to which they 
could reduce within some effective theory limit. Despite many attempts however (see e.g. [4]) an issue still under 
debate is how the spectrum of initial, pre- "Big Bang", matter fluctuations is affected by the bounce. 

To try and answer that question, a generic model for bouncing universes was studied in [f] (see also [5]) : spatially 
closed Friedmann-Lemaitre models, driven by a scalar field, whose scale factors near the bounce are expanded as a 
Taylor series in conformal time with a priori arbitrary coefficients. The conclusion of that study was that the spectrum 
of large scale matter perturbations was generically modified when going through the bounce. 

Two hypotheses underlie that conclusion. One is that the scalar field conformal gradient remains small at the 
bounce. The other provides the (required) additional information on the asymptotic behaviour of the scale factors 
(not encoded in a low order Taylor series) and is in practice, as we shall see, to approximate the scale factor in the 
bouncing region by a truncated Taylor expansion, that is by a polynomial in conformal time. 

In this Note, we first show in Section 2 that approximating the scale factor in the bouncing region by a polynomial 
amounts, when the scalar field conformal gradient is small at the bounce, to more than expanding it in a Taylor 
series. Indeed, as we shall see, the crucial fact that the effective potential for the perturbations decays at the outskirts 
of the bouncing region (and hence allows to define the transition matrix for the spectrum) is an information which 
cannot be obtained by a low order Taylor expansion but stems from the fact that the scale factor is approximated by 
a polynomial. We will however extend the class of scale factors which yield decaying potentials for the perturbations 
to other generic functions of time, which are not truncated to low order polynomials in the bouncing region. 

We then bridge a gap left open in [I] and turn in Section 3 to the characteristics of the scalar potentials required to 
reach the conclusion that the perturbation spectrum is modified as described in [I]. We find they are sharply peaked at 
the bounce and decay at its outskirts. We then study various simple scalar potentials exhibiting the same dependance 
in the scalar field at the bounce and which decay at its outskirts, and find that they also lead to modifications of the 
perturbation spectrum, but that they must be fairly fine tuned in the transition range to avoid premature recollapsc 
of the universe. 

Finally we discuss in Section 4 the hypothesis that the scalar field conformal gradient remains small at the bounce 
and argue that it can be relaxed. If, then, the scalar field conformal gradient is allowed to be of order unity at the 
bounce, then, as we shall see, either the spectrum is not modified or the universe recollapses immediately unless the 
scalar field potential remains sharply peaked. 

Section 5 summarizes our conclusions. 

II. Bouncing scale factors : Taylor expansions vs asymptotic behaviours 

This Section summarizes [I] (see also [5]), putting the spot light on the hypothesis which is made to obtain 
the desired asymptotic behaviour of the perturbation potential, to wit, the scale factors can be approximated by 
polynomials in conformal time in the bouncing region. We then extend the class of scale factors yielding similar 
perturbation potentials to other generic functions of time, which arc not truncated to low order polynomials in the 
bouncing region. 

Consider a spatially closed, homogeneous and isotropic universe with line element : ds 2 — a 2 (r])(—dr] 2 + df2§) where 
i] is the (dimensionless) conformal time, a{-q) the scale factor (with dimension of a length) and d£l\ the line element of 
a 3-dimcnsional unit sphere. If matter is just a scalar field ip with potential V(ip) the Einstein equations reduce to the 
Friedmann-Lemaitre equation : 3(H 2 + 1) = k (\p 12 + a 2 V) (k is Einstein's constant), plus cither the Klein-Gordon 
equation or, equivalently (unless <p' = 0) : 



where a prime denotes differentiation with respect to conformal time and where H = a' /a is the (dimensionless) 
Hubble parameter. In (2.1) and in the following we set Einstein's constant equal to one (so that ip and ip' become 
dimensionless) . (We also set the velocity of light to one, so that the unit of length, say, is still unspecified) . 

Consider then the perturbed metric : ds 2 = a 2 (r/)[— (I + 2^)dr/ 2 + (1 — 2^)cK7§] and the perturbed scalar field 
ip(i]) + Sp>. In Fourier space, the scalar perturbations ^ n and 5ip n are functions of time and of the Eigenvalues n 
of the Laplacian on the 3-sphere (defined as A/„ = — n(n + 2)f n , n e N and n > 2). The linearized Einstein equations 

(see, e.g. [6]) then yield two constraints ($„ = ^ n ,Sp n — 2 ^ $n , +a> "M , and the following evolution equation for 



-p' 2 =H 2 -W + 1 



(2T) 



= -2^$ 



n 




with / = 



,'2 



(2.2) 
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and k 2 = n(n + 2) — 3. According to [1], the values of k which correspond to scales of cosmological interest today lie 
in the range [60 — 6 x 10 6 ]. 

Consider now the family of scale factors indexed by the integer N : 

" ' ■ - 5(1 + « 8 = 5^. (2-3) 



a ( r l) = fl 2n — I with a = 1 , a 2 = - , a 4 = 



24 ' u 720 



As in [1] we restrict our attention to time-symmetric bounces. The value of the scale factor at the bounce is rescaled 
to be one (this choice sets the unit of length to be the radius of the universe at the bounce, see Note [7]). The 
introduction of the parameter ij allows to choose the value of a 2 at will. £ and \ are parameters, taken to be of order 
one, used in [1]. 

Near the bounce, (2.3) is the Taylor expansion to order 0{rj 2N ) of any bouncing scale factor, and yields : 



(higher order expansions will not be necessary). Clearly, r] must be greater or equal to one. As in [1] we shall not 
consider the case 770 = 1 (for which the perturbation equation (2.2) is singular at the bounce, see Note [8]). Moreover 
we shall concentrate, as in [1], on the case when the (dimensionless) scalar field conformal gradient ip' remains small 
at the bounce, that is to T <C 1 (we shall discuss this hypothesis in Section 4). 

As for the Taylor expansion of the effective potential for the perturbations W(n) it reads : 



61* 5£ f i + 23 10£ 



4ry 2 T 2 V 2 V 2 T r,^T 2 



t + Oiv 4 )- (2.5) 



With the additional assumption T <C 1, and for generic coefficients (£, \, ...) of order one, it no longer depends on \ 
and becomes, see [1] : 

W(v)^^(l + ^V 2 )+0( V 4 ) when T 0+ . (2.6) 



2T \ T 

Now, and this is one of the main points of this Note, in order to relate the pre-bounce perturbation modes to 
the post-bounce ones, one needs asymptotic regions where to define the "in" and "out" states. Such asymptotic 
behaviours cannot in general stem from a low order Taylor expansion, whose range of validity is, usually, limited : 
the Taylor expansion of W(rj) always blows up, whatever the order in r\. If, moreover, T <C 1 it generically (that 
is for (£,x 5 •■•) °f order one) blows up for 3> V / T, that is, well inside the bouncing region, see (2.6). Hence the 
assumption, made in [1] (see, e.g., equations (30) or (4-1) therein), to approximate the scale factor in the bouncing 
region by the polynomial (2.3). That extra assumption yields the following asymptotic behaviours : 



which, in practice (see Figure 1), are valid at the outskirts of the bouncing region, that is for y/T <C \rj\ < 1, see Note 
[9]. 

Concentrating mainly on the case — 1<£<0,A = 2, the authors of [1] were then able to approximate W(rj) by a 
simple rational function and show that its shape was fairly independent of £ and 770 , and hence generic. 

Therefore, with the two hypotheses T <C 1 (that is small scalar field conformal gradient at the bounce) and 
polynomial scale factor in the bouncing region, W(rf) is a deep potential well (see Figure 1), which can in effect be 
approximated by (2.6) for \r\\ < and, from (2.7), by W = for \rj\ > or, even more simply by 

a square potential of depth V m = — 5£/2T and width yj— 2Y/5£ (or, better, twice that value, see [1]). 

Today large scales correspond to modes k 2 small compared to the depth of the potential, if T is sufficiently 
small (say, T < 10~ 6 ). It is then a simple exercise to see that their spectrum is affected by the bounce. More 
precisely : the pre-bounce mode, which behaves as u n (q) = A^ e exp (ikrj) + B^ rc exp (—ikij) in the asymptotic region 
1 < K] <C — \J — 2T/5£ where W(n) ~ is the square potential approximation of (2.7), turns into the post-bounce 
mode u n (rf) = AP° st exp (ikrf) + B^° st exp {-ikt]) for 1 > 77 > y/-2T/5£_ with, (sec [1]) : 
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FIG. 1: Matter perturbations induced by polynomial scale factors. The left panel (first obtained in [1]) shows a typical deep 
potential well for the perturbations W(x) with x = rj/rjo, generated by a polynomial scale factor. The parameters are the same 
as those illustrating [1], that is : 770 = 1.01, £ = —2/5, N — 2, see equation (2.3) for definitions. The right panel exhibits a few 
perturbations modes u n {x) (generated by the potential W(x) on the left panel), and shows that their spectrum is modified by 
the bounce, according to (2.8) (see [1]) when n 2 is small compared to depth of W(x). 

where the precise value of the Const, (s TV) depends on the level of approximation made for W{rf) in the region 
\r)\ < y— 2Y/5£. (The fact that the transition matrix is not invertible is not surprising as we are in the limit where an 
incident wave is totally reflected by the deep well, see [5] for developments). The important point is the dependance 
on k, that we checked by direct numerical integration, see Figure 1, [10]. 

Other cases than — 1 < £; < 0, N = 2 can be considered. For example, we looked at the cases : N = 1 (parabolic 
scale factor), N = 3, —1 < £ < with x > or — 1 < % < : in all those examples, the shape of the potential W{rj) is 
almost independant of the higher order coefficients a^n and of the value of T (as long as it is small) , and we checked 
that the result (2.8) obtained in [1] holds (at least approximately). In fact, we found that it also holds when, say, 
N = 2, —2 < £ < — 1, that is when the universe recollapses, because, for T small enough, there is a region around 
the bounce where, first, the potential W(t]) can be approximated as above, and, second, which is large enough to 
define the asymptotic states of the perturbations u n . We also checked that scale factors polynomial in cosmic time 
t (and, hence, non polynomial in conformal time dr\ = dt/a(t)) yield the same result, that is (2.8). Finally another 
possibility is to impose that the scale factor asymptotes its quasi-de Sitter value in the range VT < \r)\ < 1, for 
example : a{rj) — 1/ cos -q + (64 i] 4 + be if ...)/ coshry. As can easily be checked, such a truncation generically has the 
same effect as choosing a polynomial scale factor : it yields a potential for the perturbations similar to Figure 1, and 
hence, once again, the result (2.8) obtained in [1] holds. 

Therefore the conclusion reached in [1], that is the modification of the perturbation spectrum through a bounce 
due to a deep potential well is quite robust, in the sense that it holds for a large class of scale factors. As we have 
stressed, the condition of validity of the conclusion are that. (1), the scalar field conformal gradient remains small at 
the bounce (T <C 1) and, (2), that the effective potential for the perturbations, W(rj), decays for \r]\ ^> yf, so that 
asymptotic regions for the in and out perturbation states can be defined. 

However, before claiming convincingly that all the examples studied "disprove a priori any general argument stating 
that the spectrum should propagate through the bounce without being modified" [1], the pertinence of the hypotheses 
underlying the claim must be probed and, hence, the characteristics of the potential V(ip) of the scalar field which 
drives these bouncing universes must be studied. 

III. The driving potentials 

The Friedmann-Lemaitre equation (which has not been used up to now) can be written as, using (2.1) : 

V(t}) = \(H' + 2H 2 + 2). (3.1) 

Near the bounce, the Taylor expansion of V(n), for any scale factor, follows from (2.3) and is, (Note [11] recalls the 
units chosen) : 

^-^-»(>-*)S + (?-f-S + S)S +o «- 
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FIG. 2: Polynomial scale factors : the driving potential. The 
polynomial bouncing scale factor (the parameters are the same 
(0 < n < 1), V(tp) can be approximated by the straight line : 
scalar field effective equation of state. 
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left panel shows a typical driving potential V(tp) which yields a 
as those illustrating [1], see Figure 1) ; in the intermediate range, 
V a {f>) = 3(1 — 0.56 <p). The right panel gives the corresponding 



Since p(j]) can be extracted from (2.4) by a simple integration, V(tp) is known. Near the bounce its Taylor expansion 
is (choosing without loss of generality ip(0) — and p' > 0, but excluding the value T = 0) : 



-wrj \ ,„ y\ (1 5£ W /2T 61* 5£ 35g 25g \ ^ 4 6 



When, as in [1], the scale factor is approximated by a polynomial, the large ip behaviour of V(ip) is : V(ip) ~ 2a 2 , 
from (3.1) ; hence V(<£>) oc rj~ AN , from (2.3), and, since tp' ~ \/2 from (2.6) : 

W ^ "° 4v + O ( -jk ) ■ (3-4) 



"2N 



p^" \(p' 



Figure 2 gives the typical shape of V(ip) when T< 1. 

The main properties of V(tp) when T « 1 arc : (1), V(0) w 3 (that is V(0) » 3/a 2 , in Planck units, see Notes 
[7] and [11]) ; (2), it is sharply peaked at ip — (as is clear from (3.3)), that is it exhibits an effective square mass 
to 2 « S which is negative and large (with respect to the scale set by the radius of the universe at the bounce ; in 
Planck units, m 2 « 4^, and may be small if ao is large enough, see Note [7]) ; (3), more striking perhaps, it turns 
out that it decreases almost linearly in the transition region before reaching its large tp behaviour (3.4). Another 
way to characterize the underlying physics is to describe the scalar field as a fluid of density g = + V((p)\ and 

pressure p = \ £p — V{p) s j and give its adiabatic index p/g as a function of g, see Figure 2. and Note [12] 

Let us now turn to simple scalar potentials which resemble these "conical hat" potentials, for example, Lorentzians 
V V: c{p) — 1+ " 2y 2 or Gaussians V v , c (ip) — ve~ c v , with v of order one and c large (we do not ask here if pre-Big- 
Bang or cyclic scenarios can provide us with such potentials). It is then an easy task to see that, when one imposes 
ip'(0) = V2T <C 1 (or, equivalently, v « 3), one gets potentials W V}C (Tf) for the perturbations which also exhibit a 
deep well and an aymptotic region, and hence for which one can also conclude that the spectrum of perturbations is 
modified by the bounce, see Figure 3. Of course, the scale factors are then no longer polynomial in time. In fact, at 
least in the examples considered, the universe soon recollapses ("soon" meaning for \rj\ « 1), because V VyC (p) does not 
decrease slowly enough in the transition region. If, therefore, one demands a model which not only yields a potential 
W{rj) for the perturbations that exhibits a deep potential well together with asymptotic regions, but which also avoids 
premature recollapse, then it seems that the potential for the scalar field must be fairly fine tuned. 

IV. Should t be small ? 

In order to avoid misunderstandings, we reestablish here the constant ao that we had previoulsy set to 1 and use 
Planck units. We also suppose ao 3> 1 so that the radius of the universe at the bounce is large compared to the 
Planck length. 

In the preceding Sections we followed [1] and imposed the dimensionless scalar field conformal gradient p' (0) = y2T, 
see (2.4), to be small. It seems to us that a more meaningful quantity is the kinetic energy density = of the 
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FIG. 3: Lorentzian scalar potentials (1). The left panel shows the scalar potential Vfjfi) = 1+< ^ 3 with v and c such that 
V(i£>) has the same small tp behaviour as in Figure 2. The central panel shows the corresponding potential for the matter 
perturbations, W(rj), which must be compared to Figure 1. (For 77 large, W(rf) diverges because the scale factor recollapses.) 
Finally, the right panel exhibits a few perturbations modes u„(x), with n 2 small compared to the depth of W(rj), and shows 
that their spectrum is modified by the bounce. 




FIG. 4: Lorentzian scalar potentials (2). The left panel shows the scalar potential V(ip) = 2 with v — 2.5 (hence T = 0.5) 

and c 2 = 10. The central panels show the corresponding scale factor. The right panel shows the potential for the matter 
perturbations, W(rf), which is similar to Figure 3. 

scalar field. At the bounce Ek\o = T/oq. Now, it turns out that, in all the cases considered (with T <C 1, see end of 
Section 2) Ek (which is even in 77) increases sharply with -q and reaches a value of order 1/a 2 , before gently decaying. 
Hence imposing Ek\o to be much smaller than its generic value in the whole bouncing region seems somewhat unatural. 

Now, an even more meaningful quantity seems to be the total energy density of the scalar field. From Friedmann- 
Lemaitre's equation, it is : 

Q (r 1 ) = ^ + V = ^(H 2 + l). (4.1) 

Its value at the bounce is g(0) — 3/a 2 ,. For a Q ^> 1, see Note [7], g(0) is small compared to the Planck density, 
whatever the value of T. Moreover, in all cases considered, g gently decreases in 0(r)~ 2N ) in the asymptotic region, 
whatever the value of T. Therefore, if one argues that what matters is the value of the total energy of the field with 
respect to the Planck energy in the bouncing region, then T need not be small. 

Now, for T (and £) of order one, the scalar field potential V(ip) is not sharply peaked at the bounce- see (3.3)-, the 
potential W(rj) for the perturbations is no longer deep compared to fc 2 -see (2.5)- and then, if there is an asymptotic 
region large enough to define the in and out states, the spectrum of the perturbations is unaffected by the bounce, as 
stated in [1]. If, on the other hand, no asymptotic region can be reached before recollapse (as is the case for "bell" -like 
potentials described by Lorentzians or Gaussians with both parameters of order one) then no definite conclusion can 
be drawn about the modification of the perturbation spectrum. 

A last possibility must be considered though : T of order one, but large (negative) £. This is a case when developping 
the scale factor in a Taylor series is not a good strategy. Such values of the parameters can however be achieved by 
simple sharply peaked potentials described by, say, Lorentzians V VtC (ip) — with 2 < v < 3, and c large. As can 

be easily checked the potential for the perturbations and the perturbations modes are then similar to those presented 
in Figure 3, the caveat being again that the decay of V(ip) must be slow enough to prevent premature recollapse, see 
Figure 4. 

V. Conclusion 

In summary, we found that the main result of [1] , that is the modification of the spectrum of matter perturbations 
through a bounce due to a deep potential well, is quite robust, provided two hypotheses are fulfilled : (1), the scalar 
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field conformal gradient remains small at the bounce (T < 1) and, (2), the effective potential for the perturbations, 
W(rj), decays for \r]\ 3> v^T, so that asymptotic regions for the in and out states can be defined. 

We also found that these hypotheses, with the additional requirement that the universe does not recollapse pre- 
maturely, yield a "conical hat" scalar field potential, with an effective squared mass which is negative and large in 
comparison with the scale set by the radius of the universe at the bounce (but which can be small in Planck units) ; 
more peculiar is the fact that the scalar field potential seems to have to decay almost linearly (before tending to zero 
in the asymptotic region) in order to prevent premature recollapse. 

Finally, we examined hypothesis (1) (T <C 1) and found that it was not necessary if one only imposed the total 
energy density of the scalar field to remain small (in Planck units) in the whole bouncing region. Now, when all 
parameters are of order one, the spectrum of the perturbations is either unaffected by the bounce, as stated in [1] 
or, if the decay of the scalar field potential V(ip) is too fast, the universe recollapses immediatly and no definite 
conclusion can be drawn on the modification of the perturbation spectrum. However a sharply peaked potential 
V((p), with 2 < V(0) < 3 (in units set by the size of the universe at the bounce), which decays slowly enough, "does 
the job" and ensures a modification of the perturbation spectrum. 

The question now is : do there exist plausible (ideally string-inspired) scalar field potentials which possess such 
properties. In this short Note we shall leave that difficult problem open. 



Acknowledgements. The author warmly thanks Jerome Martin and Patrick Peter for discussions on their paper [1], 
which led to the first version of this Note. Recent discussions on their criticisms led me to send them a second version 
of my Note, which expanded a number of points to try and make them clearer. Unfortunately that did not seem to 
have convinced the authors of [1] : see their Web Note gr-qc/0406062 It is that second version that is presented here. 
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Notes 



[7] Hence, if, say, I = 3 denotes a length that we want to compare to the Planck length Ipianck , then the dimensioned 
length corresponding to I — 3 is l p h = 3 ao Ipianck similarly if to = 3 denotes a mass, then the dimensionful mass is 
rn p h = 3mpi anc k/ao, and, if g = 3 denotes an energy density, then the dimensionful density is g p h — 3gpi anc k/aQ, etc. 
Finally, if T is a dimcnsionlcss number, then T = 3 means... T — 3. Consequently, when in the following we refer to 
a dimensionful quantity as "small" or "large", we shall mean small or large with respect to the scales set by ao- Now 
(and I thank Jerome Martin and Patrick Peter for pointing that out to me), if ao is a large number (that is if the size 
of the universe at the bounce is large with respect to the Planck length) the scales set by ao are small compared to 
the Planck scales. 

[8] In fact, in order for the perturbation equation (2.2) to make sense, ip' must be different from zero in the whole 
bouncing region. The approach developped here and in [1] therefore does not apply to the problem studied by C. 
Gordon, N. Turok, Phys. Rev. D 67 (2003) 123508 (see also N. Deruelle and A. Streich, |gr^qc/0405003|). 
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[9] There is of course nothing a priori "wrong" in approximating the scale factor by a polynomial. Indeed what 
happens is the following. Since T is supposed to be small, there are two ranges of variation for r\ : the bouncing 
region \r/\ <C 1 (in practice |r/| < 1) where the scale factor (which is supposed to be a smooth function of time) is well 
approximated by a Taylor expansion in rj (to order 4 say) , and the corresponding range of validity for the perturbation 
potential W(n), which is much smaller, |?7| < \/T. In order to obtain W(n) outside its range of validity set by the 
Taylor expansion of the scale factor (to order 4 say), more information is a priori required about the scale factor (its 
Taylor expansion to order 2N ~ 1/T say). Truncating the Taylor expansion, that is to say approximating the scale 
factor by a (4th order) polynomial is just a simple way to provide that missing information and is a perfectly valid 
thing to do if the Taylor series Lagrange remainder, that is the time derivative of the scale factor of order 5 (say) is 
small on the whole range \n\ e [0, 1]. Moreover it guarantees that the effective potential for the perturbations decays 
at the outskirts of the bouncing region, so that the transition matrix for the perturbation spectrum can be defined. 
(A particular example, studied in [1], of a scale factor which can be approximated by a truncated Taylor series of 
order 4 is a{rf) = ^fl + rf = 1 + i] 2 /2 — i] 4 /8 + 3ry 6 /16 + ... for which the Lagrange remainder is small in the whole 
bouncing region.) See however Section 3 for the implications of such assumptions on the shape of the potential V(ip). 

[10] If it happens that A^ IC — — £?P rc , then the expansion at next order in k/\/V m yields AP ost = —B^ ost = 
Const. AP rc , where the Const, is of order one. In that case then, the amplitude of the pre-bounce modes is not 
amplified and its spectrum is not modified by the bounce. 

[11] In Planck units, V(r]) in (3.2-4) must be replaced by a 2 V{rf) where ao is the radius of the universe at the 
bounce in units of the Planck length, see Note [7]. 

[12] Choosing to impose the scale factor to approach its quasi-de Sitter expression in the transition region (instead 
of being a polynomial) yields a potential V(ip) similar to Figure 2. 



